Cosmological perturbations in f(T) gravity 
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We investigate the cosmological perturbations in f(T) gravity. Examining the pure gravitational 
perturbations in the scalar sector using a diagonal vierbien, we extract the corresponding dispersion 
relation, which provides a constraint on the f(T) ansatzes that lead to a theory free of instabilities. 
Additionally, upon inclusion of the matter perturbations, we derive the fully perturbed equations 
of motion, and we study the growth of matter overdensities. We show that f(T) gravity with 
f(T) constant coincides with General Relativity, both at the background as well as at the first-order 
perturbation level. Applying our formalism to the power-law model we find that on large subhorizon 
scales (0(100 Mpc) or larger), the evolution of matter overdensity will differ from ACDM cosmology. 
Finally, examining the linear perturbations of the vector and tensor sectors, we find that (for the 
standard choice of vierbein) f(T) gravity is free of massive gravitons. 

PACS numbers: 98.80.-k, 04.50.Kd 
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I. INTRODUCTION 

Cosmological observations of the last decade [l[ indi- 
cate that the observable universe experiences an accel- 
erated expansion. Although the simplest way to explain 
this behavior is the consideration of a cosmological con- 
stant 0], the difficulties associated with the fine-tuning 
problem have led some authors to postulate more radi- 
cal alternatives. A first direction is to consider the dark 
energy paradigm, which, at least at an effective level, 
can originate from various fields, such as a canonical 
scalar field (quintessence) @, a phantom field, that is 
a scalar field with a negative sign of the kinetic term 
[H, or the combination of quintessence and phantom in 
a unified model named quintom [j|. A second direction 
is to modify gravity itself, using a function f(R) of the 
curvature scalar @, higher derivatives in the action Q, 
braneworld extensions string-inspired constructions 
@, holographic properties [l(|, UV modifications [Til ] 
etc. 

Recently, a new alternative approach appeared in the 
literature [13 - [l4l ]. It is based on the old idea of the 
"teleparallel" equivalent of General Relativity (TEGR) 
[HI, EH , which, instead of using the curvature defined via 
the Levi-Civita connection, uses the Weitzenbock connec- 
tion that has no curvature but only torsion. The dynami- 
cal objects in such a framework are the four linearly inde- 
pendent vierbeins (these are parallel vector fields which is 
what is implied by the appellations "teleparallel" or "ab- 
solute parallelism" ) . The advantage of this framework is 
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that the torsion tensor is formed solely from products of 
first derivatives of the tetrad. As described in the 
Lagrangian density, T, can then be constructed from this 
torsion tensor under the assumptions of invariance under 
general coordinate transformations, global Lorentz trans- 
formations, and the parity operation, along with requir- 
ing the Lagrangian density to be second order in the tor- 
sion tensor. However, instead of using the torsion scalar 
T the authors of [ID, EI| generalized the above formalism 
to a modified f(T) version, thus making the Lagrangian 
density a function of T, similar to the well-known exten- 
sion of f(R) Einstein- Hilbert action. 

In comparison with f(R) gravity, whose fourth-order 
equations may lead to pathologies, f(T) gravity has 
the significant advantage of possessing second-order field 
equations. This feature has led to a rapidly increasing 
interest in the literature, and apart from obtaining ac- 
celeration [l3|, EH one can reconstruct a variety of cos- 
mological evolutions [13, EH , add a scalar field [l9| , use 
observational data in order to constrain the model pa- 
rameters I20B and examine the dynamical behavior of the 
scenario |21| . 

All the previous investigations on f(T) gravity have 
focused on the background evolution. However, in order 
to reveal the full scope and physical implications of the 
theory, one must delve into the perturbative framework. 
Thus, one must first investigate the perturbations of the 
pure gravitational sector, and extract the corresponding 
dispersion relations, since such an analysis is necessary 
in order to determine the stability (or lack thereof) of 
the theory. Additionally, one should examine the com- 
plete set of gravitational plus matter perturbations, since 
they are related to the structure growth, and therefore 
their evolution could constrain f(T) gravity via compar- 
ison with the wealth of precision data from observations. 
Both of the above investigations are the main goals of 
this work. In particular, the present investigation is in- 
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terested in the first-order perturbations of a Friedmann- 
Robertson- Walker universe under f(T) gravity. 

The plan of the work is as follows. In section UU we 
present the cosmological scenario of /(T) gravity at the 
background level. In section IIIII we examine the first or- 
der perturbed equations for the scalar sector, we extract 
the corresponding dispersion relations, we derive the evo- 
lution equation for matter overdensities, and we evolve 
perturbations for a specific model as an example. In sec- 
tion IIVI we extend our analysis in the vector and tensor 
perturbations at linear order. Finally, section [V] is de- 
voted to the summary of our results. The calculations in 
the preceding sections are performed in the Newtonian 
gauge, however for completeness we include Appendix A 
with the perturbed equations for the scalar modes in the 
synchronous gauge. 

II. THE COSMOLOGICAL BACKGROUND IN 
/(T) GRAVITY 

In this section we study the cosmology of a universe 
governed by f(T) gravity. In this manuscript our no- 
tation is as follows: Greek indices fj,,i/,... run over all 
coordinate space-time 0, 1, 2, 3, lower case Latin indices 
(from the middle of the alphabet) ... run over spatial 
coordinates 1, 2, 3, capital Latin indices A, £?,... run over 
the tangent space-time 0, 1, 2, 3, and lower case Latin 
indices (from the beginning of the alphabet) a, will 
run over the tangent space spatial coordinates 1, 2, 3. 

As we stated in the Introduction, the dynamical vari- 
able of the old "teleparallel" gravity, as well as its f(T) 
extension, is the vierbein field e^^). This forms an 
orthonormal basis for the tangent space at each point 
of the manifold, that is • eg = tjabj where 
Vab = diag(l,— 1, — 1, —1). Furthermore, the vector 
can be analyzed with the use of its components e A in a 
coordinate basis, that is = e^9 M . 

In such a construction, the metric tensor is obtained 
from the dual vierbein as 

9^u(x) =r)ABe£(x)e^ (x). (1) 

Contrary to General Relativity, which uses the torsion- 
less Levi-Civita connection, in the present formalism 
ones uses the curvatureless Weitzenbock connection [22j , 

T vll = e J 4 A (9„e M A , whose torsion tensor reads 

= T vlt - = e\ (d^ - d u e*). (2) 



Moreover, the contorsion tensor, which equals the differ- 
ence between the Weitzenbock and Levi-Civita connec- 
tions, is defined as 



K »» p = __^ p _ T »» p _ Tp »uy (3) 

Finally, it proves useful to define 

S p ^ = l(K^ p + S^T^-S;T^ a ). (4) 

Note the antisymmetric relations T A = — T A „ M and 
S p ^ u = —Sp"^, as can be easily verified. Using these 
quantities one can define the so called "teleparallel La- 
grangian" as [H, H, [HJ 

L T = S p ^T^ v . ( 5 ) 

In summary, in the present formalism, all the informa- 
tion concerning the gravitational field is included in the 
torsion tensor T A , and the teleparallel Lagrangian Lt 
gives rise to the dynamical equations for the vierbein, 
which imply the Einstein equations for the metric. 

From the above discussion one can deduce that the 
teleparallel Lagrangian arises from the torsion tensor, 
similar to the way the curvature scalar arises from the 
curvature (Riemann) tensor. Thus, one can simplify the 
notation by replacing the symbol Lt by the symbol T, 
which is the torsion scalar (l4j ] . 

While in teleparallel gravity the action is constructed 
by the teleparallel Lagrangian Lt — T, the idea of f(T) 
gravity is to generalize T to a function T + f(T), which is 
similar in spirit to the generalization of the Ricci scalar R 
in the Einstein- Hilbert action to a function f(R). In par- 
ticular, the action in a universe governed by f(T) gravity 
reads: 

I = J d ' Xe t T +f( T )+ L m] , (6) 

where e = det(e^) = \J—g and L m stands for the matter 
Lagrangian. Variation of the action with respect to the 
vierbein gives the equations of motion 



e-%{eS A n^ + fin - e\T% x S p ^[l + f'(T)] + S A ^(T)f"(T) - -e A [T + f(T)} = AirGe p A T p \ (7) 

I 

where a prime denotes the derivative with respect to T and the mixed indices are used as in Sa^ = s^Sp^ . 
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Note that the tensor T p v on the right-hand side is the 
usual energy-momentum tensor, in which we added an 
overset label in order to avoid confusion with the torsion 
tensor. 

If we assume the background to be a perfect fluid, then 
the energy momentum tensor takes the form 

em 

T = V9»v - (p + pjUftUu, (8) 

where w M is the fluid four-velocity. Note that we are 
following the conventions of (25| . but with an opposite 
signature metric. In the following sections we will be in- 
terested in the perturbed energy-momentum tensor, and 
the signs of the perturbed terms are not affected due to 
the indices mixing (one upper and one lower). 

Let us now focus on cosmological scenarios in a uni- 
verse governed by f(T) gravity. Thus, throughout the 
work we consider the common choice for the form of the 
vierbien, namely 

e t = diag(l,a,a,a), (9) 

which leads to a flat Friedmann-Robertson- Walker 
(FRW) background geometry with metric 

ds 2 = dt 2 - a 2 {t) 5 rj dx l dx 1 1 (10) 

where a(t) is the scale factor. It has recently been shown 
26] that f(T) gravity does not preserve local Lorentz in- 
variance. In principle then one should study the cosmol- 
ogy resulting from a more general vierbein ansatz. In this 
paper we specialize our attention to the diagonal form, 
leaving the more general case for future investigation. 

Using this vierbein, together with the above fluid de- 
scription of matter, one sees that equations (0 lead to 
the background (Friedmann) equations 

H 2 = 8 -^f Pm - - 2f(T)H 2 (11) 

fr = AnG{p m +p m ) 

l + f'(T)-12H 2 f'(T)' [ ' 

In these expressions we have introduced the Hubble pa- 
rameter H = a/a, where a dot denotes a derivative with 
respect to coordinate time t. Moreover, p m and p m stand 
respectively for the energy density and pressure of the 
matter content of the universe, with equation-of-statc 
parameter w m = p m / p m . Finally, we have employed the 
very useful relation 

T=-6H 2 , (13) 

which straightforwardly arises from evaluation of (|5|) for 
the unperturbed vierbien ([9]). Lastly, note that Gen- 
eral Relativity is recovered by setting f(T) to a constant 
(which will play the role of a cosmological constant), as 
expected. 



III. SCALAR PERTURBATIONS IN f(T) 
GRAVITY 

One of the most decisive tests for the reliability of a 
gravitational theory is the examination of its perturba- 
tions. This investigation reveals some of the deep features 
of the theory, determining its stability and the growth of 
matter overdensities. In this section we analyze in de- 
tail the linear-order scalar perturbations of f(T) grav- 
ity, leaving the vector and tensor ones for the next sec- 
tion. In particular, we extract the full set of gravitational 
and energy-momentum tensor perturbations in subsec- 
tion HIO] and in subsection IIII Bl we examine the stabil- 
ity of the theory. Finally, in subsection IIII CI we exam- 
ine the growth of of matter overdensities. For simplicity 
we perform the calculations in the Newtonian gauge (for 
completeness, Appendix A displays the scalar perturba- 
tions in the synchronous gauge). 

A. Matter and scalar vierbein perturbations up to 
linear order 

We are interested in examining how the scalar vier- 
bein perturbations affect the equations of motion. We 
mention that while in General Relativity the fundamen- 
tal object is the metric, which is then perturbed, in f(T) 
gravity, as well as in the old "teleparallel" equivalent of 
General Relativity, the fundamental object is the vier- 
bein. Thus, the starting point is the vierbein perturba- 
tion, which will then give rise to the perturbed metric. 

Using the symbol for the perturbed vierbein and 
for the unperturbed one, the scalar perturbation writes 

< = < + (W) 

where 

= K = s l a g o=^ < = ^ (is) 

% = fyl> t^ = -S^a<p t% = -5^ #=-A0,(16) 

with indicial notation as stated at the beginning of sec- 
tion |II] As usual, we have kept terms up to first order in 
the perturbations. Moreover, unless otherwise indicated, 
subscripts zero and one will generally denote respectively 
zeroth and linear order values of quantities. 

In the above perturbations we have introduced the 
scalar modes <f> and ip, which are functions of x and t. 
These symbols, as well as the various coefficients have 
been conveniently chosen in order for the vierbein pertur- 
bation to induce a metric perturbation of the known form 
in Newtonian gauge with signature (H ), namely 

ds 2 = (1 + 2^)dt 2 - a 2 (l - 2(j})8 ij dx i dx j . (17) 

In the following, we make the simplifying assumption 
that the scalar perturbations are diagonal, a procedure 
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which is sufficient to extract qualitative results about the 
stability of the theory. Thus, the determinant becomes 



det(e^) = a 3 (l + ^-30). 



(18) 



Proceeding forward, we calculate T x ^ and to 



first order under the perturbations (|T5|) and (|16l) . The 
torsion reads 

TV = (e X A + t X A )[d^ + tt) d v (ej + )], (19) 

and thus the second is easily calculable using ©. After 
some algebra we find (indices are not summed over) 



jiO 

So 
T 



djt ; Si ij 



n, (k<t 

a 2 



T 4 0l = H — <f> 
^•fo-VO. (20) 



H 

1 

2a 



Additionally, up to first order the torsion scalar defined 
in ([5]) is found to be 

T= S P ^T\ U =T + Tx, (21) 

where 

T = -6H 2 (22) 

is the previously seen zeroth order result and 

Ti = \2H{Hi> + j>) (23) 

is the first order one. Thus, we can easily express f(T) 
and its derivatives up to first order as: 



f(T) = /(T ) + T] 



df(T) 



dT 



f(T) - 
f"(T) 



df(T) 



•2 



fo + fl 



dT 
d 2 f(T) 



T=T n 



Ti 



d 2 f(T) 



dT 2 



dT 2 
d 3 f(T) 



T=T n 



+ Ti 

T=T 1 dT 3 



T=Tr 



= fo + fl 
= fo' + f{ 



(24) 



that is the right hand sides of these equations are func- 
tions of To and linear functions of T\. 

Let us now consider the perturbations of the energy- 
momentum tensor. The perturbations are then expressed 
as 

em 

<5T ° = -Sp m (25) 

em 

5 T l = a- 2 (p m +p m )(-diSu) (26) 

em 

ST t ° = ( P m + PrnMSu) (27) 
em 

5Ti> = 5 lJ 5p m + d l d ] n s . (28) 



where 7r s is the scalar component of the anisotropic 
stress. Inserting everything in ([7]) we finally obtain 

E° ^ (1 + /^)(V 2 0) + 6(1 + / O )tf0 
+6(1 + / )#V - 3f[H 2 

-^iA = -4nG6p m , (29) 

El = (1 + f'M + (1 + f'o)Hdrt 

-\2HHf'M = -4irG{ Pm + Pm )d t Su, (30) 

El = \2B 2 d^M + ffV)/o - (1 + /o)ft«i(0 + ffVO 
= 4irG(p m + p^d^u, (31) 



£?!<5?= ^-(-3H 2 



in 

iT ) + — ( 12ii" 2 ff 



a V /a 
0(T O + f ) Ti + A 



4a 
(l + /o) 



4o 



\QH4> + QH 2 4)-'iH 2 <t) 

a 

+4> + H(2ip - 0) + i?<0] 

f" 

■12-(-24HHj> - A8ipH 2 H - 12i? 2 
a 

-12H 3 ip + 12H 2 H<p) 

4lrG i j. r \ 



(32) 



a (l + /o) 



d j S> b 3'5?(<f>-il>) 



4TrGa 2 d ] did l Sy s 



(33) 



where we have used the definition V 2 = JT did 1 and 
indices are summed over only when explicitly shown with 
the symbol. 

The above equations are perfectly general. In what 
follows, we set the scalar anisotropic stress ir s to zero 
for simplicity, a choice which precludes the possibility 
of anisotropic expansion of the Universe. However, the 
cosmological consequences of dark energy models with 
anisotropic stress have been previously considered in the 
literature [28|, and a detailed study of the effects of 
anisotropic stress in /(T) cosmology would be an inter- 
esting avenue for future research. 

The zero-anisotropic-stress assumption, which accord- 
ing to ([33]) implies <f> = ip, along with 5p m = 0, may lead 
the system of perturbation equations becoming overde- 
termined. Once the vanishing of the anisotropic stress 
is implemented, we then have four equations determin- 
ing the three remaining perturbation variables Sp m , <p 
and Su. However, in the limit /"(T) ~ 0, equa- 
tions ([3U)) and (|3"Tj) become identical, removing the over- 
determination. We therefore conclude that the require- 
ment of no anisotropic stress imposes another constraint 
on /(T) models, namely that f"(T) ~ 0. However, note 
that this requirement on /(T) might be relaxed for more 
general choices of vierbien than © . 
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B. Stability 

In the previous subsection we derived the linear equa- 
tions of motion for the perturbations. In this subsection 
we will use them to examine the stability of f(T) gravity 
and extract the dispersion relation for the pure gravita- 
tional perturbations. Since we are concerned with the 
stability of the theory, in this subsection we ignore the 
matter sector in the set of equations (|2T))) -(|3"3" )) . This is 
because if the gravitational sector itself is unstable, the 
matter content cannot cure the instability. 

Working in Fourier space we introduce the mode- 
expansion of cj> as 



d 3 k 
(2tt)I 



(34) 



Therefore, inserting this into (f23j) . we obtain the modes 
of the Fourier transformed T\ as 



T lk = \2H l <\> k + \2H<p k , 



(35) 



while ((24)) gives the modes of the Fourier transformed /i 
as 



fik = 



df(T) 



dT 



T=T 



Tikfo, 



the modes of the Fourier transformed f[ as 
d 2 f(T) 



f'lk = T lk 



dT 2 



T=T 



Tikfoi 



and the modes of the Fourier transformed /{' as 
d 3 f(T) 



flk = Tlk 



dT 3 



T=T 



Tikfi 



o i 



(36) 



(37) 



(38) 



where obviously /o, /q and ft do not depend on 0. 
Inserting equations (f34"l) -(f38 |) into (|32j) . we obtain 

4> k + T0 fe + w 2 fe = o, (39) 

where uj 2 and V are respectively given by 



uj 2 = 



^-+H-f+Hf( ) - [36H 4 - A8H 2 H) ft + UAH 4 Hft 
l + ft- 12H 2 f>> 



(40) 



AH 



r = 



1 + ft - (l2H 2 + 9H) ft + MH 2 Hft 



1 + l2H 2 ft 



(41) 



Equation (p9|) allows one to study the stability of any 
given model. In particular, a model for which cu 2 is neg- 
ative will clearly be unstable. Note that as mentioned 
before, it is sufficient to consider a scenario without any 
matter content. In such a pure gravitational case, the 
background equation (fT2|) (that is with p m = 0, p m = 0) 
leads to a constant H. Therefore, equation (14*U1) reduces 
to the simple form 



u 2 = 



l + K-umft ■ 



(42) 



One can use relation (|42[) in order to determine if a 
specific f(T) model is free of instabilities. For instance, 
we can insert the power-law ansatz f(T) = a(—T) n , 
with a = (67Jq) 1_ ™ / '(2n — 1) in the absence of mat- 
ter (Hq is the present Hubble parameter), and the ex- 
ponential ansatz f(T) = -aT (1 - e pTo/T ), with a 
1 'I — (1 — 2p)e p ] in the absence of matter, considered 
in [l4l ] , and then we can straightforwardly determine the 
allowed values for the ansatz-parameters numerically. It 
is easy to check that in both these cases the stability 
condition is satisfied for the phenomenologically relevant 



ranges of parameters, that is < n < 1 for the power-law 
model and < p < 1 for the exponential model. 



C. Growth of perturbations 

Having examined the stability of f(T) gravity, in this 
subsection we switch on the matter sector, and exam- 
ine the fluctuations about the FRW background in the 
presence of matter. This is a crucial subject in every cos- 
mological scenario, and thus it is the starting point for 
a thorough investigation of the cosmology in the f(T) 
framework. 

In order to study the growth of perturbations, we as- 
sume for simplicity a matter-only universe, that is we 
impose p rn = = Sp rn . As usual, we define the matter 
overdensity 6 as 



5p n 



(43) 
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We now write equation ([29)1 as 

3H (1 + f - l2H 2 .fJ) $ k 
+ [(3H 2 + k 2 /a 2 ) (1 + f' Q ) - 36iJ 4 /o] k 
+ AnGp5 k = 0. 



(44) 



Note that this is the relativistic version of the Poisson 
equation in f(T) gravity. In summary, equations (|39|) 
and (|44"1) can be used in order to evolve 8 for a given f(T) 
model, which will allow for a comparison to observational 
data. 

Before proceeding to the investigation of equation (14"H) , 
let us make a comment concerning the relation to General 
Relativity. As we observe, a very interesting feature that 
emerges from the above analysis, both of the pure gravi- 
tational sector of the previous subsection as well as of the 
matter inclusive sector of the present subsection, is that 
the linear perturbations in f(T) gravity reduce to those 
of General Relativity in the limit where f(T) is constant. 
For example, if f(T) — const. = — 2A, where A is the cos- 
mological constant, then one can immediately see that 
equations (|29]) - (l33|) reduce to the well-known first-order 
equations of General Relativity [25l . l27j , and furthermore 
equations and ([H]) reduce to the well-known equa- 
tions for the growth of perturbations in ACDM cosmol- 
ogy (see e.g. [2g[). 

Obtaining General Relativity as a limiting behavior at 
the level of both the background and the perturbations, 
is not a standard feature of a gravitational theory. In 
particular, any covariant modified gravitational theory 
contains extra degrees of freedom at the perturbational 
level, unless this theory can be conformally transformed 
back to the standard General Relativity. These degrees 
of freedom may have no impact at the background level 
(under certain symmetries), but the perturbative modes 
could in principle leave their imprints at first order, even 
if one imposes the General Relativity limit. This is the 
case, for example, for Hofava-Lifshitz gravity [llT ] where, 
although the theory in the Infra- Red coincides with Gen- 
eral Relativity, at the perturbative level one obtains the 
known strong coupling problem, even in the IR [30j . 

However, as we see from the analysis of the present 
work, f(T) gravity seems to coincide with General Rela- 
tivity when /(T) is constant, both at the background 
level and at the first-order perturbation level, leaving 
possible differences to arise at second order or beyond. 
Such a behavior is seen also in f(R) gravity, both in its 
metric and Palatini formulation [3l[ |32j , and it is a signif- 
icant feature of f(T) gravity and a robust self-consistency 
test. 

We now proceed to the investigation of the physical 
implications of a non-trivial /(T)-ansatz, studying the 
growth of the overdensity for a specific model. We choose 
the power-law model suggested in [l4| : 

f(T) = a(-T) n , (45) 
where a — OoHq)^ 1 (I — fl m o) /(2n — I) and H and 



il m o refer to the Hubble parameter and the matter den- 
sity parameter at present. 



k=0.1 h Mpc 



— n=0 
---n=0.2 

- - n=0.4 



to 0.4 




iog 10 (i+z) 



FIG. 1: The evolution of the matter overdensity S as a func- 
tion of the redshift z, on a scale of k — O.lh Mpc' 1 , for three 
choices of n, for the power-law model given by J^5| ). 

Our results for the growth of perturbations, arising 
from a numerical elaboration, are presented in Figures 
[T][3l In these figures we follow the growth of the mat- 
ter overdensity <5, from the time of last scattering to the 
present one, for different choices of n, for three differ- 
ent /c-scales. In Fig. 0]we depict the evolution of S for 
fixed n = 0.2, but for three different scales. As usual, we 
use the redshift z as the independent variable, defined as 
1 + z = ao/a with ciq the present scale- factor value. 

As expected, the n — case is identical to ACDM 
scenario [29| . However, as n increases we find that there 
is a suppression of growth at smaller redshifts, which can 
act as a clear distinguishing feature of these models. We 
also notice that larger scales are more strongly affected 



k=0.01 h Mpc" 



— n=0 
---n=0.2 
- - n=0.4 



«> 0.2 




iog 10 (i+z) 



FIG. 2: The evolution of the matter overdensity S as a func- 
tion of the redshift z, on a scale of k = O.Ol/i Mpc' 1 , for three 
choices of n, for the power-law model given by J^5| ). 
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,x 10 



k=0.001 h Mpc" 1 



to 5 




— n=0 
---n=0.2 

- - n=0.4 



1 2 3 

log 10 (l+z) 

FIG. 3: The evolution of the matter overdensity 8 as a func- 
tion of the redshift z, on a scale of k = O.OOl/i Mpc' 1 , for 
three choices of n, for the power-law model given by {45\ f. 



n= 0.2 



0.8 



0.6 



to 0.4 



0.2 



— k=0.001h Mpc" 1 
---k=0.01h Mpc" 1 

- - k=0.1h Mpc" 1 



12 3 

iog 10 (i+z) 

FIG. 4: The evolution of the matter overdensity S as a func- 
tion of the redshift z, for a fixed n for the power-law model 
given by [45ty , for three different scales. 



than smaller ones. Therefore, observations which target 
sub-horizon scales close to the horizon could, in principle, 
constrain these models. 

Finally, note that, as we mentioned in the end of 
section IIII A[ in order for our scenario to be free of 
over-determination (since we have imposed the zero- 
anisotropic-stress assumption which reduces the degrees 
of freedom) we must have f"(T) ~ 0. This condition in 
the case of the power-law ansatz requires n<l, which 
is what we have used for the numerical analysis. Inter- 
estingly enough, it is exactly the same condition that is 
needed in order to acquire an observationally compatible 
dark-energy and Newton constant phenomenology at the 
background level [l4| (for the power law ansatz n ~ 1 
also gives f"(T) ~ but we do not consider this case 
due to the aforementioned phenomenological reason). 



IV. VECTOR AND TENSOR PERTURBATIONS 
IN f(T) GRAVITY 

In the previous section we focused on the scalar pertur- 
bations of f(T) gravity, since they are sufficient to reveal 
the basic features of the theory, allowing for a discussion 
of the growth of matter overdensities. For completeness, 
in this section we extend our analysis in order to include 
the vector and tensor sectors of the theory in the absence 
of matter. 

The general perturbed vierbein at linear order reads 



as 



SKl + ^+aiGi + diF)^ 

ck5J(1 -</>)+ a {hi a + d l d a B + diC a + d a C t ) <5* 

#(1-V0 

~ [<*£(! + 4) + (h\ + d l d a B + d l C a + d a C l ) <5f ] 



(46) 



In these expressions, apart from the scalar modes <f> and 
if) of the previous section, we have introduced the trans- 
verse vector modes Gi and Ci, the transverse traceless 
tensor mode hf, and the scalar modes F and B, the di- 
vergence of which will also contribute to the vector sector. 
Similarly to the simple scalar case, the coefficients on the 
above expressions have been chosen in order for this vier- 
bein perturbation to give rise to a perturbed metric of the 
familiar form: 



.9oo = 1 + 2-0 

9m = a [diF + Gi] 

gij = -a 2 [(1 - 2^)6, 



didjB 



diCj] 



(47) 



Let us make a comment here concerning the number of 
degrees of freedom of the perturbed theory. As may be 
deduced straightaway, , KZ V and S ptlu are spacetime 
tensors under an infinitesimal coordinate transformation 
of the form 



(48) 



This implies that the torsion scalar T is a generally co- 
variant scalar, and thus actions of the form of (JBJ) are 
generally covariant as well as invariant under (|48p . As 
a result, for our choice of vierbien ©, the number of 
degrees of freedom (DOF) is identical to General Rela- 
tivity. In particular, in 3 + 1 spacetime dimensions, the 
metric, being symmetric, has 10 independent DOF. This 
is reflected in (|4"6"|) , which comprises 

• 4 scalar DOF ip, 4>, F and B 

• 4 vector DOF, 2 associated with each of the diver- 
genceless vectors Gi and C, 
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• 2 tensor DOF associated with the transverse, trace- 
less and symmetric tensor hij 

However, not all of these DOF are independent as there 
exist 3 + 1 DOF associated with the coordinate transfor- 
mation (|3S]l (the temporal part of is the scalar e , and 
its spatial part can be decomposed into the gradient of a 
scalar plus a divergenceless vector: di€ S + ej ' , leading to 
a total of 2 scalar and 2 vector DOF). Subtracting these, 
we are left with a total of 6 DOF: 2 scalar, 2 vector, and 
2 tensor, just as in the case of General Relativity. 

We can therefore work in the Newtonian gauge, setting 
F and B to zero. This is easily understood since under 
(|48p. the gauge transformation of B is — (2e s )/a 2 , while 
that of F is (l/a)(-e - e s + 2He s ). Therefore, e s can 
be chosen in order to give rise to B — 0, and similarly an 
accompanying choice of eo will lead to F = 0. 

From the above analysis, it is clear that there are no 
extra modes in f(T) theories, which often show up in 
theories with less symmetry than General Relativity. For 
example as discussed in [33| , in the case of Hof ava grav- 
ity, which is invariant under spatial diffeomorphisms and 
space-independent time re-parametrizations but not un- 
der space-dependent time re-parametrizations, one can 
no longer choose the longitudinal gauge, and thus an ex- 
tra scalar DOF remains. However, note that extra DOF 
can still arise from more general choices of vierbien than 
©. The possible existence and phenomenology of these 
extra modes is beyond the scope of this work and war- 
rants further study. 

Additionally, we choose a gauge where the vector mode 
Ci vanishes (through an appropriate choice of ej). As 
usual, the vector modes are transverse, while the tensor 
mode is transverse and traceless, namely 

d,C l = d t G l = ; dih ij = 5 l] h l3 = 0. (49) 

Finally, we easily deduce the following relation between 
the tensor perturbations in the vielbein and inverse viel- 
bein: 

h p= i a=1 = -h a=1 »=\ (50) 

Using the above relations, the perturbed torsion tensor 
([H?f becomes 

T\ v = d^-d^Sl + a^Gv-dvGv) 
T% = drf + dih^Si-tyhfSi, (51) 

where for notational compactness we have introduced a 
Go part of Gi, which is zero. This torsion tensor inserted 



into Q leads to 

S 0t = ^+^EW+- G ' 
a 1 2a z ' a 

3 



S ij = 


4a 3 y 3 


- did) 


= 


-H + cj) + 


2H4, + ±k a 6l 


Si i0 = 




- 8^) 


Si ij = 








+-Gj + 
a 


a 

2a 3 



However, the torsion scalar is unaffected by the vector 
and tensor modes, and thus it remains 

T = T + Ti = -6H 2 + 12H 2 tP + 12H<j), (53) 

and likewise the determinant e is still given by 

e = a 3 (l + r/» - 30). (54) 

We now have all the necessary machinery in order to 
extract the equations of motion for the vector and ten- 
sor sectors. Following the steps of the previous section, 
we can similarly decompose the energy-momentum ten- 
sor into its vector and tensor components and ignore the 
vector and tensor anisotropic stresses. We finally obtain 

[1 + J'(T)]V 2 G J = (55) 

for the vector mode. Since the quantity in square brack- 
ets is zero only for the unphysical model /(T) = — T (for 
which the action © does not describe the gravitational 
sector anymore), we can eliminate it, resulting in 

V 2 Gj = 0. (56) 

Therefore, the vector modes in f(T) gravity decay as 
1/a 2 , that is similar to the General Relativity case. 
For the tensor mode, we obtain 

(57) 

Finally, similar to the scalar case, we can Taylor-expand 
the derivatives of f(T) using (|24|) . and Fourier-expand 
the vector and tensor modes, in order to obtain the corre- 
sponding dispersion relations. Moreover, we can split the 
tensor sector into left-handed and right-handed polariza- 
tions. However, such a detailed analysis of the gravita- 
tional wave spectrum lies beyond the scope of this work. 
Here we retain only the simple forms (|56j) and (|57|) , since 
they are adequate in order to reveal the basic features of 
the behavior. 

Concerning the tensor equation (|57l) , although there is 
a new friction term, there are no new mass terms, which 
is a behavior similar to the scalar case of the previous 
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section. Therefore, we can safely conclude that, in gen- 
eral, f(T) theories do not introduce massive gravitons; 
thus when f(T) tends to a constant we do not obtain the 
typical problems of massive gravity, which is a significant 
advantage of f(T) gravity. Additionally, note that simi- 
lar to the scalar case, in the limit where f(T) tends to a 
constant we do recover the behavior of General Relativ- 
ity at linear order, which is a self-consistency test of the 
construction. Lastly, from the equations of motion for 
scalar, vector and tensor perturbations presented above, 
it is clear that these three classes of perturbations decou- 
ple from one another in f(T) gravity, just as they do in 
case of General Relativity. 



and deserves further investigation. 
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Appendix A: Scalar Gravitational and matter 
perturbations of f(T) cosmology in the synchronous 
gauge 



V. CONCLUSIONS 

In this work we investigated the recently developed 
f(T) gravity, going beyond the simple background level. 
f(T) gravity is the extension of the "teleparallel" equiva- 
lent of General Relativity, which uses the zero curvature 
Weitzenbock connection instead of the torsionless Levi- 
Civita connection, in the same lines as f(R) gravity is 
the extension of standard General Relativity. In partic- 
ular, we analyzed the first order perturbations of f(T) 
gravity. 

Examining the scalar perturbations of f(T) gravity in 
the Newtonian gauge, we derived the perturbed equa- 
tions of motion and extracted the corresponding disper- 
sion relation. Therefore, in constructing a realistic f(T)- 
gravitational scenario, one should use an f(T) ansatz 
that leads to non- negative to 2 in (j42|), in order to ob- 
tain a theory free of instabilities. Moreover, we showed 
that for the assumption of no scalar anisotropic stress to 
be consistent, one needs the constraint f"(T) ~ 0. 

Additionally, we found that f(T) gravity with f(T) 
set to a constant coincides with General Relativity, not 
only at the level of the background but also for the first- 
order perturbations. This is a significant advantage of the 
theory as compared to other modified gravity paradigms. 

Furthermore, as an example of an application of our 
formalism, we followed the growth of perturbations in a 
specific f(T) model, namely the power-law ansatz pro- 
posed in [14(. For this model we found that on large 
subhorizon scales (0(100 Mpc) or larger), the evolution 
of the matter overdensity differs from ACDM. Therefore, 
future precise observational data on these scales could be 
used to constrain or rule out such models. 

Finally, we investigated the vector and tensor per- 
turbations at linear order, extracting the corresponding 
equations of motion for the vector and tensor modes. 
We showed that f(T) gravity does not introduce mas- 
sive gravitons, which is a significant advantage. Lastly, 
we verified again, as in the scalar sector, that in the limit 
where f(T) tends to a constant the theory tends to Gen- 
eral Relativity, both at the background as well as at the 
linear perturbation level. 

Clearly f(T) cosmology presents a very rich behavior 



Let us analyze the scalar gravitational and matter per- 
turbations in the synchronous gauge. The perturbed 
vierbein read 



e l = e % 4- f 



(Al) 



with 



p° — t)° ■ p a — fr a n ■ p m — f)^ ■ — — 



, = n • r = ,v 



a 



*°=0 ; t; = 5 b ^(5 a b A-a 2 d b d a B) 
tf=0 ; tZ = -6^{6 b a A-a 2 d b d a B) 



(A2) 



These symbols, as well as the various coefficients have 
been conveniently chosen in order for the vierbein per- 
turbation to induce a metric perturbation of the known 
form in synchronous gauge, namely 

ds 2 = dt 2 - a 2 [5ij(l + A) + didjB] dx l dx } . (A3) 

Inserting these into the perturbed torsion tensor (|19[) . 
and then to (U]) we obtain 



T\ v = SI \ H+^- l -d t {a 2 d l d b B)5 b v 

T'ji = \d 3 A 

S ° l = \d l A 

S?i = \&A 
4 

Si 0j = ~^d t (a 2 d^B) 

Si 0i = - +\d t (a 2 J29 j S b j d i 6 i b B) 



H + - J + -d t {a 2 \7 2 B - a 2 9 i 5 i B).(A4) 
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Furthermore, up to second order the torsion scalar de- Finally, concerning the first-order matter perturbations, 

fined in ([5]) is found to be we use the results (|25 p -(|28 |) . 

, Inserting these results into ([7]) we extract the per- 

T = -6H 2 - -PAdiA - 6HA + 2Hd t (a 2 V 2 B). (A5) turbed equations of motion as: 



E° 



{ 1±^V 2 A 3H 2 f[(T) - 3HA(1 + Q + Hd t {a 2 V 2 B){l + f Q ) Tl±MH 



= AttGSTq , 



K = 



E 



El 



-d l A(l + f Q ) + Hd l A(l + f Q ) - 6HHf[l(T)d*A = AttGS T a \ 

^ di 8iA + ^m di 8i 
la a 

AttG „ em 



i J2 ! d ! :.\i),S:.\ - 6HA + 2Hd t (a 2 V 2 B) 



(A6) 
(A7) 



(A8) 



.^±fld t [a 2 d t (a 2 d 3 Sl^B)} + { l±I^d. } 5l^A + + / ° ) d t {a? d^d* B) { ±±J^d t {Ha^d 3 8lB) 



fl!(T) ■ ■ a 4-77(7 em 

.$i*±J. HHdt tf d 8>PB) + ea&HfgWPdjSlB - -c^5(T + /„) = ST 3% 

a o a 



(A9) 



3H 2 + H 



(l + M^^dti-Ha 2 ) (- 



ZA 



a 2 V 2 B 



-HA - ^ + ^d t (a 2 V 2 B - a 2 d l d l B) + y {a 2 V 2 B - a 2 d l d l B) 



dt^a 2 

— (V 2 A -did 1 A) 
4a 

HA+^d t [a 2 (W 2 B + d t d'B)^ ' r{{T)X2H 2 R 



fS(T) 



H 



d t (J2 d l Ad l A) - 6H 2 A - 6H 2 HA 



HHd t {a 2 V 2 B) + 2H 2 d 2 {a 2 V 2 B) + 6H 2 Ha 2 d l d l B + 3HHd t (a 2 d i d i B) 

m±A) + (To + /o) (^_£ ft3 . B 

em 1 n 

4ttG5 T 1 1 - 47rGp m —A + inGp^d^B, 
la I 



(A10) 



r 



where we have used the definition V 2 = Y^i®i® l > anc ^ the S symbol, 
indices are summed over only when explicitly shown with 
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